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Recent measurements suggest our universe has a substantial dark energy component, which is 
usually interpreted in terms of a cosmological constant. Here we examine how much the form of this 
dark energy can be modified while still retaining an acceptable fit to the high redshift supernova 
data. We first consider changes in the dark energy equation of state and then explore a model in 
which the dark energy is interpreted as a fluid with a bulk viscosity. 

PACS numbers: 98.80 -k, 98.80.Es 


1. INTRODUCTION 


Publication of the analysis of type la supernovae red shift data by the High-z Supernova Search Team Q 
and the Supernova Cosmology Project 0] provided the first indication that the universe is accelerating. The 
data favor an interpretation in which most of the energy in the universe is some form of dark energy capable 
of providing negative pressure. The usual candidate for the dark energy is a cosmological constant A, which, 
as the name suggests, gives a time independent dark energy density. Both groups analyze their data assuming 
contributions from ordinary matter and a cosmological constant. They present their results as contours the 
the plane of the parameters UIm and Da, defined below, which are a measure of the relative contributions of 
matter and dark energy QQ From these plots (e.g., Fig.^and TableP), it is clear that values of UIm and Da 
are most probably those of an accelerating cosmology. 

While this picture is attractive because of its simplicity, one can ask to what extent the data can distinguish 
between alternative hypotheses for the behavior of the dark energy. Investigations into this question are generally 
framed in terms of determining the equation of state of the dark component (A), that is, determining w = px / px , 
wher^x is the pressure of the dark component and px is its energy density, as a function of the scale parameter 
R{t) |3,3- If i® possible to quantify the effect of any particular w by fitting the supernova data and comparing 
the resulting value of with the one for w = —1, which corresponds to the case of a cosmological constant. 

In the next Section, we review briefly the strategy for extracting the cosmological parameters Dm and Da from 
the effective magnitude data for type la supernovae, and use a straightforward minimization to reproduce 
the results of Refs. Q and Q. In Section 0 we use the values of y^ for w = — 1 to explore several modifications 
to the equation of state including a search for the minimum y^ for the case w = constant. Changes in the 
pressure-densiW relation of this t ype or similar modifications, including models using a Chapl ygin equation of 
state 3, B 0, lE 0, llE Il4 flR or generalized Cardassian expansion [3 llE: have been 

considered before. Our discussion is included here as a prologue to Section^] where we present a model in which 
the dark component is a fluid with a bulk viscosity. It is shown that this model provides an equally good fit to 
the supernova data, predicts an accelerating universe and has mass density fluctuations which grow. All of this 
occurs at the expense of some entropy production. Finally, we conclude with some comments on the current 
supernova data’s capacity to discriminate between models of the dark energy. 


2. FITTING THE SUPERNOVAE DATA WITH A COSMOLOGICAL CONSTANT 


If R{t) is the scale parameter in the Robertson-Walker metric and k its curvature constant, the Friedmann 
equations with a cosmological constant are |^. 
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where G is Newton’s constant, p is the density of matter, p its pressure and a dot denotes time differentiation. 
Conservation of energy and momentum gives the additional relation 

d{pR^) =-pd{R^). (3) 


Letting = R^/R^, and denoting the present time by the subscript 0, Eq. can be divided by iLg, the 
square of the Hubble constant, to give 
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in general. Assuming the pressure of ordinary matter is negligible, Eq. © gives 
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and, in terms of the scaled variables y = R/Rq and r = Ho{t — to), Eq. @ becomes 
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From Eq. 0 , the time history of the universe can be obtained as 
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Owing to the cubic behavior of the denominator of Eq. (jS]) , the Hm — Ha plane is divided into different categories 
of histories depending on whether a zero of the cubic occurs for ?/ < 0, 0 < ?/ < 1, or y > 1. These are, 
respectively, those universes which continue to increase in size from y = Q, those which never achieve ?/ = 0 and 
those which start at y = 0, increase to a maximum size and then collapse. The boundaries of these regions can 
be seen, for example, in Fig.^along with the line corresponding to flat cosmologies, Ha = 1 — Hm- Also, at the 
present time, Eq. © can be written 
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and hence the universe is accelerating if Ha > HAf/2. 

The extraction of the cosmological parameters Hm and Ha is achieved by relating the apparent magnitudes 
of m{z) of the type la supernovae and their absolute magnitudes M to their luminosity distance dL{z) as 

m(2;) - M = 51ogio ((iL(z)) + 25 , (10) 

where c^l is measured in megaparsecs. The luminosity distance of a supernova at red shift z can be expressed as 
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with c denoting the velocity of light. In Eq. cu, sin[ ] is used for ilk > 0, sinh[ ] is used for ilk < 0, and the 
unmodified square bracket is used for ilk = 0. 

To obtain a fit to the data, we write the apparent magnitude as 


m{z,ilM,ilA,M) = ^logiQiVL) + M , (12) 

where = c~^HodL and Ai = M + 5logiQ{cHQ^) + 25 is an additive constant. The parameters Hm, and 
Ai are then determined by minimizing 
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with (Ti denoting the error in mexp(zi)- 

Our fit to the data presented in Ref. is shown in Fig.^ To avoid repetition, we restrict our fits to 
this data. We have obtained similar results for the data in Ref. P|. When computing points on the = 
— Xmin contours, we choose the minimum value of Ai at each point (Om,11a)- This is easily done because 
the minimization condition for Ai can be solved explicitly in terms of Hm and Ha- The values of Ay^ plotted 
(2.30,4.61,6.18) would correspond to 67.3%, 90.0% and 95.5% confidence contours if the errors were Gaussian. 
The contours are not precisely elliptical indicating non-Gaussian behavior, and this is addressed in the detailed 
analysis of Refs. Q and |^. Fortunately, our more naive analysis yields very similar Ay^ contours, which we can 
use to assess the quality of the fits resulting from different assumptions about the dark component. 

For the most part, the fit shown in Fig.^ lies above the dashed line separating accelerating from non¬ 
accelerating cosmologies. While the best fit is not a flat cosmology, the y^ values in Tabled show that requiring 
Hm + Ha = 1 does not dramatically alter the quality of the fit. 


3. VARYING THE EQUATION OF STATE 


To generalize beyond the case of a constant dark energy, we can write Eq. o as 
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and determine the dark energy density px using the energy conservation equation 

dR 


dpx = -3 {px +px) 

together with the equation of state 

Px = w{R) Px . 

For models of this type, the generalization of Eq. 0 is 
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When w{y) is a constant, say w = —a, a > 0, the boundaries in the Hm — Hx plane between regions of 
expanding universes and universes which eventually collapse or have no big bang are obtained by determining 
where the minimum of the function 


Hm — {Hm + Hx — l)y + Hxy^°‘ 
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vanishes. This leads to the equation 

Ox (3q-1)(3°-1) (0^_l) N3a 

Om (3a)3“ Om / 

from which Ox can be found as a function of Om- When o = 1, Eq. gives the familiar results for the case 
of a cosmological constant |2^. For these models, the generalization of Eq. © is 
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At the present time, this gives the acceleration parameter 
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and hence the universe is accelerating if 
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To gain some sense of what varying the value of a does to the quality of the fits, we performed a analysis 
of the data in Ref. |2| for the case a = —2/3 and for the case where a is varied along with the other parameters 
to obtain the best fit. The results are shown in Figs.|^ andJP and the best fits are given in Table in These 
examples show that the quality of the fits to the data of Ref. for the cases considered - w = — 1, —2/3 , —0.73 
- is basically the same. Even when the value of w is allowed to vary, the minimum in is not lowered in any 
significant way. One may argue that values of u> > —2/3 are less preferred since, as suggested by Fig.|51 their y^ 
contours encroach into the region of cosmologies which have no big bang. 

The constant equations of state discussed above {w = —o) all result power law dark energy densities of the 
form 

= y3(a-l) ^ (24) 

pxo 

and all lead to fits with comparable y^. To determine how changing the functional dependence of px on y affects 
the quality of the fits, we examined equations of state of the form w = —by, b > 0. These equations of state lead 
to dark energy densities with an exponential dependence of the form 


px _ 
Pxo 


(25) 


In order to make direct comparisons with the constant w cases, we fit the data of Ref. 0 using b = 21^ and 
6 = 1 . These fits are shown in Figs.^landO The numerical details are given in Table lTm Again, the y^ of these 
fits indicates that data are not too sensitive to the ^-dependence of w. However, the location and size of the 
contours suggest that the w = —y choice is the more desirable. 


The lifetime of the universe depends on the choice of w and can be used to discriminate among various models. 
Fig.El shows the Hm — ^x plane for four different choices of w along with lines of constant lifetime. Generally, 
the y^ contours lie between 11.9 and 19.0 Gyr. 


4. DARK ENERGY AS A VISCOUS FLUID 


As an alternative to simply modifying the equation of state of the dark component, in this section we consider 
the consequences of assuming the dark energy consists of a fluid having a bulk viscosity (/. If a bulk viscosity 
term is introduced into the energy-momentum tensor of an ideal fluid, one effect is to replace Eq. (j21l with |2l| 


d{px R^) = -p*xdiR^) = 


Px-3C^U{R^). 


(26) 
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In view of Eq. d, the expression for R/R depends on the dark energy density px in a nonlinear fashion, which 
means that, unlike the case of the multiplicative equations considered previously, Eq. d is not separable. In 
general, it must be solved numerically. 

If we assume that the pressure px of the dark component is negligible, the differential equation determining 
Px is 


dpx = -3 


Px 
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since Ti. is the quantity needed to determine the luminosity distance d^. This differential equation is 
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The second equality follows from Eq. d- Since the luminosity distance is given by 
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where 7d(z) is obtained by replacing y by 1/(1 + z), we write the differential equation in terms of z, and solve 


dTp{z) Ti? {z) — ATL{z) 

dz (1 + z) 


f^fc(l + z). 


(33) 


From the form of Eq. d. the expression for the apparent magnitude m{z) depends on the parameters A and 
flk, which we can vary to obtain a fit to the supernova data. We perform the integration in Eq. d numerically, 
determining 7i(z') using a four point Runge-Kutta algorithm and the boundary condition 7i(0) = 1 to first solve 
the differential equation. Given di^(z), for the data point z^, we then minimize ns in Eqs. 11211 and lldll . The 
resulting contours are shown in Fig.[3and the best fits are listed in Table ITVI In Fig.[7| the lines corresponding 
to various values of Hm indicate the boundary of the region below which, for that Hm, the dark energy density 
Px remains positive. This is determined by the condition 


VLx— = 7p{z) - (1 + zf^M - (1 + zfnk > 0 , (34) 

Pxo 

which is analyzed by choosing value of SIm , fixing A, taking flk to be large and negative and solving Eq. I■‘l■lll for 
1 < z < 1000. If condition Eq. 18411 is satisfied for all z, ilk is increased and the process is repeated until, for a 
sufficiently large value of ilk, Eq. d is violated for some z. This determines a point (A, fife) on the curve for 
the given Hm and the region above this curve is not allowed. 
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The range of lifetimes for this case, obtained from the relation 


toHo 



1 

y'^iy) ’ 


(35) 


is shown in Fig.|Sl Here, unlike the cases of Secs. |2] and |21 the lifetime cannot be much less than 14 Gyr. The 
lines of constant lifetime terminate on the dark black line in this figure and Fig.[7| For and A above this line, 
the lifetime integral does not converge. In each of these figures, the dotted line separates the accelerating from 
the non-accelerating cosmologies, with the region above the line corresponding to accelerating cosmologies. This 
can be seen by noting that analog of Eq. m is 
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Hence, accelerating cosmologies satisfy ilk > —3H-I- 1. 

When p and p are related by an equation of state of the form p = w p, the cosmology evolves at constant 
entropy. In the present case, there is entropy production. This can be seen by noting that 
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where T{R) is the temperature of the dark energy. On using Eq. (12611 . we have 
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where H^y) is obtained by solving Eq. ll^ . Setting 

T-i'^iy) = y~^u{y) 

and using Eq. m to eliminate C, the change in the entropy density at present is given by 
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where pcrit = 3Hq/8ttG is the critical density and u{y) satisfies 
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To complete the calculation of ASq, the functional form of T{y) is needed. The choice of T{y) must be consistent 
with the integrability of the entropy, which, for T taken to be a function of V and px , requires |25| 
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Using Eq. (gnj) and Eq. (gnii, it is possible to obtain the relation 
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Assuming px = 0, Eq. gni implies that T(V, px) has the general form 

T{V,px)=T{Vpx)= J d/3A(/3)(Upx)^ 


(46) 


(47) 


and, for the purpose of assessing the entropy growth, we examine a single term of the form 

T = C{Vpx)^, 

where C is a constant. In this case, Eq. (HU) becomes 

dT _ %H dR _ 3A^/yy/u{y)dy _ d[u{y) - Um - y^k] 

T px R u{jj)-^lM-y^k [u{y) - ^Im - y^k] 

where Eqs. (P|l . (03) and lEll have been used. Integrating Eq. 03 gives 
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which is consistent with Eq. 14811 . Using Eq. 15011 . Eq. 14811 is also integrable and the entropy density at present is 
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As with the luminosity distance above, we evaluate u(0) numerically using a simple four-point algorithm to 
integrate Eq. (Ell . The ratio of the present value of the entropy density produced by the viscous medium to the 
entropy density of the microwave background is plotted in Eig.|n|for several values of (3 and T(l) taken to be the 
temperature of the microwave background, = 2.73K. The numerical results scale as T^/T{1). Clearly, a side 
effect of the viscosity being adequate to produce an accelerating cosmology is the production of large amounts 
of entropy. 

Although we are assuming that the matter and dark energy densities constitute non-interacting components 
of the total energy density, it is nevertheless necessary to check the compatibility of this assumption with the 
existence of matter density fluctuations which grow at large z (or small y). Eor an inviscid fluid in the Newtonian 
limit, a matter density fluctuation 5 satisfies the equation 
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6 + 2-d-4TrGpMS = 0, (52) 
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and the influence of the viscous dark energy enters through the R/R term. Using Eqs. 03 (03, the 

differential equation Eq. 03, and assuming that pM oc (1 -|- z)^, the equation for S can be written 

(1 + z)H^{z) 6" + i {n^{z) - 3AH{z) - Ufc(l + z)") - ^Um(1 + z)^S = 0, 
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where the ' denotes a derivative with respect to z. To solve this equation for large z, we note that Eqs. ( 021 ) and 


(I44II imply 

lim H'^iy) ^ y ^ (■u(O) + y) . 
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This gives for large z 

l~L^{^z) (1 + 2:)^'u(0) + (1 + z)^ . 

(55) 

Keeping only the leading terms for large z, Eq. 15311 becomes 



2{l + zf5" + {l + z)5’ = 

u(0) 

(56) 


and a solution of the form (5 = (1 + z)“ gives 

= ( 57 ) 

Hence, the solution which grows as z decreases is 

In general, a; < 1 because the condition that px be positive, Eq. m, implies, in view of Eq. ( 1 ^ . that u( 0 ) > Hm- 
The curves of constant exponent x for several values of Hm are shown in Fig.^| The regions below the x = 1 
curves have matter density fluctuations which grow, the trend being toward slower growth for larger negative 
values of Slfc. 


5. CONCLUSIONS 


When analyzed in terms of a universe consisting of matter and dark energy in the form of a cosmological 
constant, the supernova data of Refs. Hi favor an accelerating universe. The best fits in this case do not 
correspond to a flat cosmology, but a flat solution is well within the likelihood contours with a relatively close 
to the minimum value of y^. Our ability to reproduce the results of Refs. 0 is indicated in TableD 

As a means of assessing the sensitivity of the currently available supernova data to the form of the equation 
of state p = w{y)p, we repeated the analysis using equations of state of the form 

w{y) = —a — by a, 5 > 0. (59) 

For 6 = 0 and a number of different values of a, with 0 < a < 1, we find best fits with y^ very nearly identical 
to those obtained for a cosmological constant. This includes a determination of the best value of a by treating 
a as one of the search parameters. In these cases, the sum VLm + tends to be a bit larger than + Ha, 
and the flat cosmology value of Hm is somewhat smaller than the corresponding value obtain for a cosmological 
constant. When a = 0 and b is varied, the quality of the fits in terms of y^ is again virtually unchanged. The 
case w = —y is very similar to w = —1 in all respects, as can be seen by comparing Figs.^andlS] This tendency 
of the minimum value of y^ to be the same as long as the equation of state describes a dark energy with negative 
pressure is another indication that the precise form of w is ill determined by the current supernova data [23 ■ 

To examine whether the introduction of a negative pressure component by means of a multiplicative equation 
of state was essential, we explored the notion that the dark energy is a fluid with bulk viscosity. The bulk 
viscosity provides a negative pressure and, again, a fit to the supernova data gives a y^ that is indistinguishable 
from the cosmological constant case (compare Tables IIVI and IJ) . We have shown that such a dark energy can 
provide a fit which is consistent with an accelerating flat cosmology having reasonable values for the lifetime 
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and and density fluctuations with the appropriate behavior. There remains a question of whether the use 
of a simple linear bulk viscosity term Cv/y is valid throughout the range of j/(t)’s probed by the supernova data. 
There are known non-linear effects, usually associated with viscosity-driven inflation models, which could modify 
our results [^1^1^ . We have not attempted to include them in this investigation. The model presented here 
is different from those which sugges t that viscosity effects associated with cold dark matter could be the source 
of the observed acceleration . Viscosity effects have also been discussed in models with matter creation 

M 

This approach to the inclusion of bulk viscosity is highly testable in the sense that future measurements could 
rule it out. As mentioned above, the lifetime cannot be much less than 14 Gyr. Further, if the universe is flat, 
or nearly so, as seems likely, then we require to be rather small. A close examination of Fig. linear 0^ = 0 
reveals that the present supernova data rule out = 0.3 at the 90% level. On the other hand, improving the 
supernova contours is unlikely to rule out this model, even for a flat universe, as long as Om < 0.1 is allowed. 
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Supernova Cosmology Project 



General 

Flat 

General 

Flat 

Qm 

0.79 

0.281^8 

0.73 

0.281^8 

Qa 

1.41 

1 — Qm 

1.32 

1 — Qm 

M 

23.9 

23.9 




56.9 

57.7 

56.0 



TABLE I: The first pair of entries corresponds to the fit obtained using Eq. m and the second pair are the results from 

Ref.|2. 



w = 

General 

-2/3 

Flat 

Qm 

0.54 


o 
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00 

1 + 

o o 
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00 o 

Qa 

2.33 


1 — Qm 
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23.9 


24.0 


56.81 


58.2 


w = 

-0.73 



General 


Flat 

Qm 

0.63 


n 1 o+o.ii 

U.iO_o.io 

Qa 

2.10 


1 — Qm 

M 

23.9 


24.0 

X^ 

56.80 


58.1 


TABLE If: Fits to the data of Ref. 0 for two constant values of w obtained using Eq. m are shown. 



w = —2yj3 
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M 

23.9 

24.0 

X^ 

56.9 

58.0 


w = -y 



General 

Flat 

Qm 

0.55 

0.21±0.08 

Qa 

1.39 

1 — Qm 

M 

23.9 

23.9 


56.9 

57.6 


TABLE III: Fits to the data of Ref. 0 for the cases w = —2j//3 and w = —y obtained using Eq. II .‘-{II are shown. 
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Bulk Viscosity 


General 

Flat 


-1.36 

0.0 

A 

1.47 

0.66 

M 

23.9 

24.0 


56.9 

58.0 


TABLE IV: Fits to the data of Ref. |2 for the case of a dark component consisting of an ideal fluid with a bulk viscosity 
are shown. 





FIG. 1: The contours for the fit to the 54 supernovae analyzed in Ref. are shown. The contours correspond to the 68.3%, 
90% and 95.4% confidence levels. The cross indicates the location of the minimum and the dot indicates the location of the minimum 
assuming a fiat cosmology. The dashed line separates accelerating from non-accelerating cosmologies. 
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FIG. 2: The Ax^ contours for the fit to the 54 supernovae analyzed in Ref. Qj, assuming a constant equation of state of the form 
w = —2/3, are shown. The contours correspond to the 68.3%, 90% and 95.4% confidence levels. The cross indicates the location of 
the minimum and the dot indicates the location of the minimum assuming a fiat cosmology. The dashed line separates accelerating 
from non-accelerating cosmologies. 



Om 


FIG. 3: The Ax^ contours for the best fit to the 54 supernovae analyzed in Ref. P|, assuming a constant equation of state of the 
form w = —a, are shown. The contours correspond to the 68.3%, 90% and 95.4% confidence levels. The cross indicates the location 
of the minimum and the dot indicates the location of the minimum assuming a fiat cosmology. The dashed line separates accelerating 
from non-accelerating cosmologies. 
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FIG. 4: The Ax^ contours for the best fit to the 54 supernovae analyzed in Ref. [3, assuming an equation of state of the form 
w = —2yl2>^ are shown. The contours correspond to the 68.3%, 90% and 95.4% confidence levels. The cross indicates the location of 
the minimum and the dot indicates the location of the minimum assuming a flat cosmology. The dashed line separates accelerating 
from non-accelerating cosmologies. 
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FIG. 5: The Ax^ contours for the best fit to the 54 supernovae analyzed in Ref. [3, assuming an equation of state of the form 
w = —y, are shown. The contours correspond to the 68.3%, 90% and 95.4% confidence levels. The cross indicates the location of 
the minimum and the dot indicates the location of the minimum assuming a flat cosmology. The dashed line separates accelerating 
from non-accelerating cosmologies. 
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FIG. 6: The contours and (dashed) lines of constant lifetime are shown. Starting from the lower right and moving to the upper 
left of each panel, the dashed green lines correspond to 8.0, 9.5, 11.9, 14.3 and 19.0 Gyr. In computing the lifetimes, we used Ho=63 
km/sec/Mpc. In the w = —2yl3 panel, the 19.0 Gyr line is not visible on the portion of the plane shown. 
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FIG. 7: The Ax^ contours for the best fit to the 54 supernovae analyzed in Ref. |3, assuming a viscous medium, are shown. The 
contours correspond to the 68.3%, 90% and 95.4% confidence levels. The cross indicates the location of the minimum th® 

dot indicates the location of the minimum x^ assuming a flat cosmology. The dashed lines are the boundaries of the region below 
which for the given value of th® dark energy density px in positive. The region above the dotted line corresponds to accelerating 
cosmologies. 



A 

FIG. 8: The Ax^ contours for the best fit to the 54 supernovae analyzed in Ref. |3, assuming a viscous medium, are shown. The 
contours correspond to the 68.3%, 90% and 95.4% confidence levels. The cross indicates the location of the minimum x^ and the 
dot indicates the location of the minimum x^ assuming a fiat cosmology. The dashed lines are the lines of constant lifetime. The 
region above the dotted line corresponds to accelerating cosmologies. 
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FIG. 9: The ratio of the present entropy associated with the viscous medium to the present entropy of the microwave background 
is plotted for the 68% contour corresponding to the best fit to the 54 supernovae analyzed in Ref. |3 a,nd = 0.1. Results 

for the indicated values of the parameter 0 are shown. The cross indicates the location of the minimum and the dot indicates the 
location of the minimum assuming a flat cosmology. The dashed portions of the curves correspond to the = 0.1 boundary 
in Fig.EI 





FIG. 10: Curves of constant values of the exponent x which controls the growth of the matter density fluctuations 5^ 5 = 1/(1+ 2 )^, 
are shown for several values of 











